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I. INTRODUCTION 

Gravitationally bound quantum states of ultracold neutrons have been discovered in the last decade jTHS]. These 
quantum states are formed if ultracold neutrons with very little energy are brought on top of a horizontal mirror that 
has essentially infinitely high potential for the neutrons. The early experiments were sensitive to the shape of the 
wave function of these quantum states. Higher precision is achievable in gravitational resonance spectroscopy, that is, 
if the differences in energy of the lowest gravitational quantum states are detected [B]. The interaction that couples 
different quantum states can be a vibration of the bottom mirror, or an oscillating magnetic field. The former is 
planned in the QuBounce project, and has shown first results CUE]. The latter is planned for the GRANIT project 
[MEI]. The ultimate goal of the GRANIT collaboration is to perform gravitational resonance spectroscopy on ultra¬ 
cold neutrons trapped in quantum states m- The length of the observation time would allow a very high precision 
in the measurement of energy differences between those quantum states. 

The motivation for these measurements comes from the sensitivity of the gravitationally bound quantum states to 
extra short-range interactions (See [13] and references therein) that might be spin-dependent, or from the possibility 
to detect Chameleons [T4]. Furthermore, the study of gravitationally bound quantum states allows for a test of the 
weak equivalence principle that is probably the only one that uses bound quantum states, albeit its sensitivity cannot 
compete with tests using free-falling cold atoms m or macroscopic bodies m- 

The paper is organized as follows: In section 2, we introduce spectroscopy of quantum mechanical bound states, 
and describe the setups that have been proposed to detect them. In section 3, we introduce the Stern-Gerlach shift 
that gives the largest shift in the setups that use a spatially varying magnetic field. In section 4, we introduce the 
spectator state shift, which is present in all measurement schemes. In section 5, we discuss the interference shift that 
is washed out for all measurement schemes that use a time-dependent perturbation. 


II. SPECTROSCOPY OF QUANTUM-MECHANICAL BOUND STATES 

The one-dimensional quantum states of ultracold neutrons in the gravitational field and above a perfect reflecting 
horizontal mirror at z = 0 are given by ip m (z) = A m Ai((z — z m )/zo ) (for z > 0), where Ai(z) is the Airy function, 

•I/O 

~o = (fi 2 /2 rn^g) ' ~ 5.87 //m is the characteristic length scale of the problem, and z m is the m th solution of 
Ai(—z m /zo) = 0. The normalization factor can be written as A m = 1/ |Ai 7 (— z m /zo)\. Energy eigenvalues of the 
lowest quantum states are very small, E m = m n gz m , and are on the order of peV. The height of the lowest quantum 
states is of the order of tens of micrometers. 

A generic setup to perform spectroscopy consists of the following steps: 

(1.) Preparation of an initial quantum state. 

(2.) Transition to the final quantum state using a periodic potential. The transition is only efficient if the frequency 
of the periodic potential matches the energy difference between initial and final state. 

(3.) State analysis and detection. 

There have been several proposals to achieve gravitational resonance spectroscopy in flow-through mode, that is, 
using ultracold neutrons that traverse a setup: 

(A) Magnetic transitions in DC mode: A generic setup is shown in Fig. |T] The downward-going step in the 
bottom mirror was proposed first in ref. HZI as a tool to depopulate the ground state. A set of parallel wires 
with an oscillating current pattern generates a magnetic field that rotates multiple times; a possible choice is 
shown in the picture. In addition, there is a vertical holding field. This resulting rotating magnetic field gradient 
induces transitions between quantum states if the rotation frequency (as seen by the neutron) matches the energy 
difference of the quantum states. A ground state filter made from the bottom mirror and a scatterer with a rough 
bottom surface that leaves an open slit height of about 2 s i; t ~ 25 /jm accepts ground state neutrons and rejects 
neutrons in higher quantum states. Earlier experiments m demonstrate the operation of this filter. The 
resonance condition in the transition region can only be fulfilled for neutrons with a certain velocity in forward 
direction. After a short free-fall region, neutrons fall a height which is given by their velocity component in 
forward direction. The neutron detector at the end of the free-fall region is position sensitive in the vertical 
direction, and a given vertical coordinate corresponds to a given forward velocity. This setup has been proposed 
earlier in mm- 

(B) Magnetic transitions in AC mode: The setup is similar to the one in DC-mode, but here, the current in 
the wires is oscillating in time. The vertical holding field is replaced by a very small holding field along the 
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FIG. 1. Sketch of the flow-through setup. Neutrons enter from the left. They go through the state preparation region (1.), 
transition region (2.), state analysis region (3.) and detection (4.). The choice of parameters used in the simulations in this 
paper is a current pattern Iq, Iq/\/2 ,0, —7o/\/2, —7o, —Iq/V2, 0,... that repeats 24 times, and a periodicity of A = 1 cm. 


wire direction, which serves only to avoid spin flip transitions at the zero crossings of the current. The neutron 
spin is rotating with the magnetic field, and the temporal oscillation of the field magnitude leads to transition 
if the frequency matches an energy difference between quantum states. The neutron velocity matters only as 
it influences the time the neutron spends in the transition region, and therefore the neutron detector does not 
need to be position-sensitive. This setup has been described in more detail by Pignol et al. m- 

(C) Mechanically induced transitions: An alternate way to produce an oscillating potential that induces quan¬ 
tum state transitions is to vibrate the bottom mirror. A proposal has been made in m for such an experiment, 
with the additional feature that the long transition region is split into two short ones with a long bottom mirror 
in between, similar to Ramsey spectroscopy in atomic physics. Before the transition region, a state filter selects 
only the ground state. Behind the transition region, another state filter selects the ground state again, and a 
detector behind the setup would count only the neutrons for which no transition has occurred. Results have 
been reported from a more condensed version of this setup by Jenke et al. im where the three regions have 
not been separated; we are not discussing this complication here. 

The Hamiltonian for an UCN above a horizontal mirror at z = 0 in the transition region equations is 

H | ip(t)) = H 0 | ip(t)) + V (t) | ip{t)) = ih ^ | ip(t)) (1) 

with 

H 2 d 2 

H ° = -^ a ? +m " 92 ' (2) 

V(t) is the periodic potential applied to induce the transitions. For magnetically induced transitions, V(t) = —jl ■ B. 

The eigenvectors of the solution to the Schroedinger equation with Hamiltonian H 0 are | m). With transition potential, 

the solutions can still be written as a superposition of these eigenvectors, but with coefficients that are time-dependent. 

Neutrons that start out to be in the m th vertical quantum state will undergo transitions to other vertical states. 

We use the ansatz 


I ip(t)} = ' nEmt H • ( 3 ) 

m= 1 

We define = ( Ei — E m )/h. The coefficients a m (t ) have to fulfill the coupled equations 

, . OO 

■ (4) 

H 2=1 

We note that V ( t ) needs not only to oscillate in time, but also to contain some dependence on z; otherwise it would 
not couple different quantum states. An oscillating uniform magnetic field would not induce transitions between 
quantum states. 
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III. STERN-GERLACH SHIFT 


In this section we discuss a frequency shift due to the magnetic field configuration in the setups (A) and (B) described 
above. The oscillating magnetic field is given by a combination of a magnetic holding field B a and a magnetic field 
that is produced by a system of parallel wires shown in Fig. [T] and which oscillates in the rest frame of the neutron. 
In appendix [A] we show how to compute the magnetic field as a function of x and z. 

In both setups (A) and (B), the neutron spin is following the magnetic field direction adiabatically. For setup (A), 
this has been shown in appendix |b| For setup (B), this has been shown in (18]. Hence, we know that —/7 n • B = 
+s±g n • \B\, where s± = +1 for "spin-up" neutrons, and s± = —1 for "spin-down"-neutrons (the magnetic moment 
of the neutron points opposite to its spin). Through this paper, we will use g n = |/t n |, which makes /z n a positive 
number. We can disregard the spin state motion, as the spin stays aligned with the magnetic field, and does not affect 
the computation other than through the sign in s±. We can give a simplified Schroedinger equation: 


H \i/j(t)) = H 0 | ip(t)) + s±n n \B\ 


IV’W) = i^ t H(t)) 


( 5 ) 


A magnetic field component that oscillates, but does not depend on z, does not introduce transitions, as the matrix 
elements (m\s±fj, n \B\\l) vanish for m^l . Therefore, we need to expand \B\ at least up to the linear term in z: 


B(z, t) 


B(z = 0 ,t) 
B(z = 0,t) 


d 

+ Z d~z 


B(z, t) 


B x (z = 0,t)^ + B z (z = 0,t)^ 


dB z 


+ Z- 


B{z = 0,f) 


\j( b x sin Lot'j + ^ B 0 + B z cos Lot^J 


(^—B x sincjf^ (—/9a; sinut) + ^ B 0 + B z cos uot^J (3 Z cos oot 
^ B x sin ootj + ^ B 0 + B z cos cot^ 


We use this in the magnetic part of the Hamiltonian, s±g n \B\, and expand in Fourier components. We get: 

s±[i n \B\ — s±fi n (a 0 + ot\ cos cot + cos 2 ut + .. .) +s±/r n (/3q + Pi cos Lot + fa cos 2 Lot + ...) z + ... 


• ( 6 ) 


( 7 ) 


|£(0,f)| 


l\ B M\ 


2=0 


Out of these, we neglect terms that oscillate with frequency 2 uo and above. We disregard the term s±[i n ao'. This term 
is constant in space and time, and it just gives a constant contribution to the energy. Furthermore, we disregard the 
term s±fi n ai cos Lot for reasons that will become clear later. The term s±/j, n faz cos Lot depends on z, and oscillates 
with lo. This term is the primarily responsible one for quantum state transitions. In setup (A), in the high field limit 
discussed in Refs. Ison, it reduces to s±fi n faz cos oot. For a lower holding field Bo, it is reduced, as shown in eqs. 
|A21| In our discussion, we use B 0 = 1.5 mT, and the field values given in appendix [A] 

The term s±fi n PoZ is not oscillating, and can be absorbed into a redefined unperturbed Hamiltonian Hq.± in ^ 
through a redefinition of g: 


H 0 .± = - 


with g ± = g + s± 


H 2 d 2 
2 m n dx 2 
Mn , 


m n g±z 

~ 3 (1 ±0.125) 


( 8 ) 

( 9 ) 


Consequences of the redefined gravitational acceleration are new basis vectors (|m±)) to new energies E rn ± = 
Em.(g±/g) 2 ^ 3 ■ Numbers given here and in the remainder of this section are for setup (A), where the transition 
of interest is 3 —^ 1. The resonance frequency := (± 3 ,± — E\,±)/h is given by 


w 13,± — w 13 1 + s± 


Mn 

rn n g 


fa 


2/3 


( 10 ) 


We note that there is a shift in the resonance frequency whose sign and relative magnitude does not depend on the 
quantum states involved in the transition: 


Awi3 = Wl3.± — w 13 


+238 rad/s 
—249 rad/s 


for "spin-up" 
for "spin-down" 


( 11 ) 
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For an unpolarized neutron beam, in setup (A) and (B) one will encounter two distinct peaks for each transition 
to ++ l, corresponding to the resonance frequencies for "spin-up" and "spin-down" neutrons, respectively. A precise 
measurement of the energies of the gravitationally bound quantum states from a single peak can be obtained if the 
magnetic field is known well enough to determine the size of the frequency shift. However, better than relying on the 
precision of a magnetic field map, is to take the following average that is not sensitive to the frequency shift: 


W 13 


U! 


3/2 

13,+ 


■ U) 


3/2 ' 
13,- 


2/3 


( 12 ) 


In this way, the precision of the measurement is not limited by this effect. 

The cause of the frequency shift is the combination of magnetic holding field and rotating magnetic field that causes 
a magnetic force on the neutron spin which does not cancel when averaged over the rotation period of the magnetic 
field. We want to call this frequency shift the Stern-Gerlach shift. In first order, it is proportional to the magnetic 
moment and a magnetic field gradient, and it changes sign with the spin orientation, as seen from 

_ 2 Mn a 1 ( Mil a ^ 4 ( L Ln Q \ i 

- — q s ±-A) - n -A) ) + s± -I -/?o 1 + ■ ■ • 

wi 3 3 m n g 9 \m n g ) 3 \3m n g J 

_ ± MnAA ( Hn/3 X B X \ ± -3+ 2g, n /3 z B 2 x B z + 2/i n /3 X B X B 2 

3 m n gB 0 y6 m n gB 0 I 24 m n gB 0 3 

±1 ( „ ^ | +... =±0.0784-0.0015 ±0.0044 ±0.0001 -... . (13) 

6 y3m n gB Q J 

If we use Hq,± in the Schrodinger equation ([5|, we arrive at 

a 

H\il>(t)) = H 0 ± \ip(t)) + s ± fi n /3 1 z cosut \ip(t)) = ih— \ijj{t)) . (14) 

at 

We find its solution with the following ansatz: 

oo 

IV’(A) = a m( t ) e ~* Bm,±t \ m± ) • ( 15 ) 

m—1 

The coefficients a^(t) need to fulfill the coupled equations 


dt 


-T- 

1=1 


a ± e -iu, ml+ t 


±t (m=b | s±fj, n PiZ | Z=b) cos cot 


(16) 


We retain only the first and third coefficients, as transitions between states are only possible close to a resonance, 
and the resonances are distinct. We will validate this approximation later. 


e“ 13,±t (3± | s±iJ, n /3iz \ 1±) cos cot — *-af (3± | s±g n /3iz | 3±) cos ut 

UjL lb lb 

f e _I “ 13,±t (1± I s±g, n f3iz I 3±) cos ut — f af (1± I s±g, n [l\z 11±) cos ut 
dtn n 


We call the terms that couple a basis state with itself the self-coupling 


terms. The coefficients are: 


Hi± := \ (1± | s±n n fi\z | 1±) 
a 

fi 3 ± := \ (3± | s±g n /3 1 z \ 3±) 
li 


+394 rad/s for "spin-up" 
—428 rad/s for "spin-down" 

+930 rad/s for "spin-up" 

— 1012 rad/s for "spin-down" 


(17) 


(18) 

(19) 


For this computation, we have used (to ± |z|m±) = (to|^|to) ( g±/g ) 1/,s , and eq. 
terms is a Rabi-Frequency, modified by the Stern-Gerlach shift: 


(C31. 


The coefficient of the first 


n 


_ i 

3i,± — r 
n 


(1+ | s±n n fd\z | 3+) 


—50rad/s for "spin-up" 

54 rad/s for "spin-down" 


( 20 ) 
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Note that the Rabi frequency limit for high magnetic holding field, = s±/r n /3 z (11 z | 3) = =F63.7rad/s, differs by 
the use of unmodified eigenfunctions, and by the reduction of fd\ to j3 z . 

To be able to remove the self-coupling terms, we substitute exp(if2 m= i- sin ojt/oj). The new fulfill 


da. 


—— = — e 1 ^ 13 '* 4 cos wt ■ e* 


dt 

daf 

dt 


= —*fi 3 i,±af e lu>13 ' ±t cos cot ■ e l 


( 21 ) 


We note that any additional term in the interaction that gives an equal contribution to and exactly vanishes. 
That allowed us to disregard s±fj, n a\ cos cot in eq. Q: This term adds to the self-coupling, but its contribution is 
independent of the state and does change nothing but a global phase. Our next step is to average over fast-oscillating 
terms, that is, the terms that oscillate with ui or faster. This is called the "rotating wave approximation" in atomic 
physics: Above, we have two kinds of these terms. First, we replace the last factor on the right side of the last 
equation system with its average value: 


41 — 1 — - 140 - 1 - 

- — sin ut 


1 - 


(^3± — ^1±) 2 

4w 2 


( 22 ) 


In Ref. [20], the quality of this approximation has been studied. The result is that to first order in the size of the 
neglected term, there is no frequency shift of the observed resonance. Their consequence is only that the solutions 
show small oscillations with frequency ui , which we also see in our numerical simulation of eqs. (C 2 ). 


Close to resonance (to ~ oj 3 i 1 ±), the right side of eq. (22 1 evaluates to about 99%. This is too small to be observed. 
Therefore, in the rotating wave approximation, the self-coupling terms can be neglected. 

The second use of the rotating wave approximation is that we write cos cot = (exp (iuit) + exp(— iivt))/2 and retain 
only the slowly varying component in the differential equation system. We arrive at: 


da 3 
dt 
da^ 
dt 


-^3i,±e- <( ‘ ,, -“ ls ' ±)t af 

-^31 


(23) 


The effect of the fast rotating components has already been studied in Ref. [2T. They lead to the so-called Bloch- 
Siegert shift in the observed resonance frequency. The Bloch-Siegert shift is too small to be observed in our experiment. 


The general solution to eq. (231 is 


if (t) 

if(t) 


U 33 U 31 
U 13 Un 


(0) 

if(0) 


(24) 


The coefficients u m i are given through 


U33 = 


COS ( J (w - Wi3,± ) 2 + ^ 3 1,± f ) + * 


w - Wi3,± 


--—= sin ( \]{U - Wl 3,±) 2 + ^ 31 ,± 0 

(w — Wl 3i ±) + fi§! ± 




U31 = 


u 11 = 


n 


31, ± 


1 — &13,±) + ^ 31 ,± 


: Sill 


(^J (w - Wl3,±)“ + ^31,± 9 ^ 


- 5(4; —Wi 3 ,±)t 


e 2 


cos ('y / '(w-a;i 3 i± ) 2 -|-n§ 1 _ ± | N ) -i “ = sin (yj(w - aq 3;± ) 2 + A') 

V 2J ^/(w-a; 1 3, ± ) 2 + nl li± ^ 




U13 = 


n 


31, ± 


: Sill 


(w — Wl 3i ±) + 


^ _ w 13,±)“ + ^ 31 ,± 7^ 


(w 13,±)* 


e ^2 


(25) 
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The evolution of the ground state population is described by a modified Rabi formula, where the modification is 
the use of a spin-dependent resonance frequency and a spin-dependent Rabi frequency that are both 

different from their value at high magnetic holding fields. In particular, for a neutron that starts in state 3 (that is, 
a 3 (°) = !,af (0) = 0), we get: 


p 3 ^ 1 = U(t)\ 2 = 


sin 2 ^ \J(lo - to i3 ,±) 2 + |^ 


1 + 


( ± y 

V ^31, ± ) 


(26) 


For comparison with our simulations, we apply an additional correction. We note that the next term in the power 
series for magnetic field magnitude, eq. © , is 


z 2 d 2 

Yd4 


B(z,t ) 


$xzB x 

4 B 0 4R 0 


(27) 


In first order perturbation theory, the terms shown lead to an additional energy shift of each state. Using eq. (C4|, 
we get 


(m± | D ± z 2 | m±) = ^ z 2 m {g±/g ) 2/3 

Px*B x 

4B 0 + 4 B 0 


(28) 


with D± = s±fi n 


66T/m 2 +6oT/m 2 


We correct the resonance frequency for that energy shift 

w 'i3.± = w i3,± + ^ (4 - z i) (9±/9) 


8 Mn-P± f 7 2 __ r 2 ) (ri: /n )- 2/3 _ / +9 rad/s for "spin-up" 

15 h 


— 11 rad/s for "spin-down" 


(29) 



FIG. 2. Ground state population for a neutron that is initially in state three and with "spin-down", for different frequencies. 
A resonance is found for u> ~ 2660 rad/s. The line denoted "modified Rabi" is developed in this section. 


Fig. [2]shows the result of a numerical solution of eqs. ( |C2| ) for neutrons that are initially "spin-down", and in state 
three. The simulation includes the quantum states one to five, and both spin states. If the rotation frequency w of the 
magnetic field is chosen to be close to a resonance frequency, population is transfered between the quantum states at 
resonance. The rotation frequency is chosen to introduce transitions between states 1 -o- 3, and the figure shows the 
appearance (and later disappearance) of the ground state population |ai-(-| 2 + |an.| 2 . For the optimum frequency, the 



























ground state population oscillated between 0 and nearly 100%. Further away from the resonance, the ground state 
population stays low. In addition, Fig. 2] shows the expectation from eq. (26) (dashed line). The analytic function 
reproduces the numerical simulation well. In the high field limit, the maximum ground state population is achieved 
after t = 7 r/ST 31 = 0.049s on resonance, that is, for a resonance frequency of W 13 = 2910rad/s. We observe that the 
population change at resonance is slower for our magnetic held than it would be for a high magnetic holding held, 
consistent with the expectation of a Stern-Gerlach shift that we just developed. 



mag. field rotation frequency co [rad/s] 


FIG. 3. Maximum transition probability as a function of the rotation frequency of the magnetic field u. The full (black) bullets 
on the left are obtained in a numerical simulation of eqs. ( |C2| ) for "spin-down", and the open (red) bullets are for "spin-up". 
The lines through the bullets give the expectation from the modified Rabi formula. For comparison, we show in the center 
(blue) the expectation for high magnetic field, that is, without Stern-Gerlach shift. 


In Fig. [3] we show the maximum transition probability into state 1 for neutrons that are initially in state 3, and 
"spin-up" (black) or "spin-down" (red). The maximum is taken over a time longer than an oscillation period. The 
observed Stern-Gerlach shift confirms our analytical estimate in eq. ( |29| ) in sign and magnitude to about 0.1%. An 
interesting outcome of this discussion is that the through-going neutrons are polarized, with opposite spin states for 
the two horizontal velocities classes that belong to the two peaks in Fig. [3] 

In addition, Fig. [3]shows the expectation for high magnetic holding fields (Bq > 30mT) that is discussed in |5Mll|: 
in this limit the Stern-Gerlach shift is negligible, and the dynamics of the system does not depend on the spin state. 


IV. SPECTATOR STATE SHIFT 

Another type of frequency shift is what we call the spectator state shift. The issue is that the rotating magnetic 
field couples not only the two states that are in resonance. The coupling to other energy states leads to a resonance 
frequency shift that is present in any measurement scheme. The spectator state shift is is formally similar to the AC 
stark shift in atomic physics. 

In the following, we will use the formalism of quasienergies, based on Floquet’s theorem. We can write solutions of 
Eq. (Q, in position space, in the form 


^{z,t) = ^2 °™e ' hErrit fm(z,t ) . (30) 

m =1 

Here, E m are quasienergies. The theorem guarantees that f m (z,t) is periodic in time, with the same period as the 
periodic potential, We can expand this function in a Fourier series: 

OO OO 

m= 1 k ——00 


(31) 
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Using V(z,t) = cos ut ■ U(z ), the quasienergy harmonics 4>m( z ) have to satisfy 

(E m + fc^)^fy)=ff 0 ^fy) + ^C/fy)(^- 1 fy)+^ +1 fy)) . (32) 

This is an equation system for each positive integer m, and for each integer k. The amplitudes a m are found from 
the initial condition. For the most simple case, that is, for neutrons that begin in state |no), the amplitudes can be 
written as 


= J2 { ( t )k m\ n o) • 


k =—oo 

Thus, the transition amplitude T no _ >ni from initial state no to final state ni after time t turns out to be 


OO OO 


Tn 0 ^n i — { n l 

m—1 k,k '=—oo 


^k' \ jrEmt—ik'cut / jji 


{4> k m I no) 


(33) 


(34) 


It is instructive to see how the Rabi formula can be obtained with restricting the sum to the two quasi-harmonics 
no < n\. Such an approximation is justified for weak (|f2 nortl | <C |w„ oril |), but resonant coupling of levels no and ni, 
that is with a detuning S := w — ( E ni — E no )/h which is small (|5| <C |w noni |). Eq. (32 I is approximated by 


E n J° no (z) = Ho<t>° no (z) + \u{z)4>\ 0 (z) , 
(- En 0 + Hoj) 0i o (z) = Hofa(z) + l -U{z)4>l 0 (z) . 


(35) 


We are looking for a zeroth order solution for which 4>n 0 ( z ) ~ Co ( z \ no) and <f>h 0 {z) ~ C\ (z \ ni). With this ansatz, 
and neglecting higher order terms, this equation system can be turned into the algebraic equations 


{^Enp E no ^j Co — 2 ^nin 0 Cl 

i^Enp T EjJ E n -^j C\ — -^nonxCo 


(36) 


As defined previously, htt nino = (no|U|ni). This system of equations has two solutions for quasienergy E no and 
coefficients Cq and (fy: 


E np _E n Q (5 1 

h h 22 

T7,n 


\J + |u„ ir 


c ± 

and —= 


C, 


o S± 


\J f 5 2 + |U„ 


(37) 

(38) 


We determine the coefficients through the normalization condition f |</>° o (z) + e lLJt (f>h 0 {z)\~ dz = 1: 

S± ' “ 


Co = 


^Js 2 + \n r 


^(s±^s 2 + + \n nino \ 2 V ' + 


^ !1± 


(39) 


With Eq. (331, we arrive at 

|T(z,f)> - C+e~^ ( C+ |no) + e~ iut C+ |m» + Cfye-^-V (C 0 ~ |n 0 ) + e^Cf |m» 


COS 


\/s 2 + \n ni , 

^nino 

s 2 + \n nino \ 2 


X - * 


sjs 2 + |n. 


: sin 


\A 2 + |U„ 


nin 0 | 




\ n 0 ) 


sin + |U„ irio | 2 0 e *( + 2 )‘|nfy 


(40) 
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We find the complete time-dependent solution to be the same as given in Eq. (24 ). If we neglect spectator states, 


resonance is achieved if the frequency of the perturbation potential oj coincides with the energy difference u>„ 


We 


will now show that taking into account additional harmonics in the quasienergy equation system (321 results in a shift 
of the resonance line in second order in Q TM nn /u ) nnn ,. For this, we need the first order correction to the wave function. 


We can compute this correction using Eq. (35). Together with the zeroth order, we obtain: 


4>n 0 ( z ) ~ °o i z I n 0 ) + Y C o,k {z | k) 
k^n 0 

4> l n 0 {z) « Cl {z I n x ) + Y C'l.fe ( z I k ) 


with 


^ _ 1 (k\U\nx) 

L-0,fc — -— 1^1 


fa k 


“ E„. 0 — E k 


2 E n o — E k 


Cx 


Cx ,k = 77 


(k\U\n 0 ) 


^ Enn + fau — E k 


-Co 


fa 0. 


nok 


2 E n i — E k 


-Cn 


Again we simplify equation system (321, but we will keep two more harmonics: 


( E no - Hu,) (k~l(z) = H 0 ^~l (z) + \u(z)<t>° no (z) , 

Eu 0 4>n 0 {z) = Ho4>° no {z) + ^U(z)</)-°{z) + (z^Yz) 
(E no +fau) c^ 0 (z) = H 0 cj)\ o (z) + ^U(z)<j>l 0 (z) + 1U(z)4>l o {z ) 
( E no + 2 Hw) <j>l o (z) = H 0 <j>l o {z) + '-^U{z)<j>n Q ( z ) 


(41) 

(42) 


(43) 


(44) 


We can turn this equation system into two coupled equations using the Green’s function approach. Greene’s operator 
G(e) is 




(45) 


Using Greene’s operator, we write the two new harmonics as 


We obtain 


<f>n 0 ( z ) = 2 G ( En o - fau)U{z')4>l 0 {z') 
and 4>l 0 (z) = \c(E no + 2fa u)U{z')^ no {z') . 


E n oCW = (Ho + -U(z)G(E no - faw)U(z') + -U{z)^(z) 


and ( E na + Hu,) = (h 0 + \u{z)G{E no + 2 Hu,)U{z')) <,(*) + \u(z)ct> 0 no (z) 


(46) 


(47) 


This equation system is similar to (351, with the only difference that we have an additional effective potential of the 
type V = (1/4 )UGU. This equation system includes the second order corrections due to spectator states. We need 
to find the solution with all terms of second order in fL 


l/w. 


Using Eq. (142]) , we obtain 


n Q n l / UJ ' n 0 n l • 


^ E Uq EnQ ^11^ Cq — 2^ n i n o Cl ^ ^ ^ ^/cno ^1,/c , 


k^n i 


{^Eno + EjJ E ni V22^ C\ — ^ ^noni^Q ^ ^ ^ 


k^no 


(48) 
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We have used the definitions 


Vn = - ( n 0 \UG(E no - huj)U\n 0 ) 


and V 2 2 = 7 (ni\UG(E no + 2huj)U\ni) 


A 2-^ 


in 


kriQ 


k =1 


h 2 


2-E-no — -E'm — E k 

l^fcnj 2 


4 ^ 2E ni — E no — Ek 


(49) 

(50) 


Substituting equation system (43) in equation system (48), we obtain 


^E nQ E nQ Vn Ull'j Cq — O ^nino^l 


Co = ^Slr. 


h , 


{Eno + hoj — E ni — V 22 — E/ 22 ) Cl — — Ori 0 m Cq 


(51) 


We introduced the notation 


4 kjt m En 0 +tujj-E k 


\(no\U\k)\< 


H v 

I 




l^fcnp | 

E„. — E k 


U 22 = T E 

4 fc^n 0 ^"0 _ 


\(ni\U\k)\ 


H v ■> 

4 


k^no 


Evjr. E U 


The solution for the quasienergies becomes: 


h 


E; 


np 

h 


6 + (Vn -V 22 + Uu -U 22 )/ft J_1 [77777 7 77t TT \ /tW , 10 

---± -y (0 + (Vir - V22 + u11 - u 2 2 ) /h) + 


(52) 


(53) 


We can compute the quasienergies E[ o , and the low-frequency part of the wave function | 4 '( 2 ,t)) J by substituting <5 
with S + (Vii — V 22 + Un — U 22 ) /ft in Eqs. (37) and (401. The maximum probability for disappearance from initial 


state no is achieved for a resonance frequency that is larger than the energy difference of the two states in resonance 
by <5max = — (V n — V 22 + Un ~ U 22 ) /h: 


A - K 

Cmax — , 


\^kni | 


4 l fe _ 1 2 E ni — E no — Ek 


E 

ky£n 0 


00 


\^kno | 


| ^4 kn 


l^fcrai r 

E no ~” Ek Jyy] 2 E no — E ni — E k E ni — E k 


(54) 


We note that the maximum probability for appearance into final state n\ can be slightly different. If we had only the 
two states no and n\, we would find (i max = (7i/4) |fl Jl0 „ 1 |" /(E ni — E n0 ), which, in this approximation, corresponds 
to the Bloch-Siegert Shift [2T. 

We now discuss the spectator state shift in setup (A) and setup (B) in the high magnetic holding field limit, that is, 
without taking into account the Stern-Gerlach shift introduced in the previous section. It can be re-introduced simply 
by treating each spin state separately; for each spin state, one would take the effective values for gravitational quantum 
state energies and Rabi frequencies. Table [I] compares the high-field limit of the spectator state shift for different 
transitions from Eq. (54) with the frequency shift from a numerical simulation using many coupled states. For each 
transitions, the strength of the periodic potential was chosen so that the time for maximum transition ir/Q nino = 0.1 s. 
The sums are converging only slowly, which means, that even quantum states with an energy considerably higher 
than E ni + fno contribute to the result. The results of the analytical expression (541 and the numerical simulation 
are in reasonable agreement. The agreement could be improved by using an iterative procedure, as we have used 
unperturbed values of energy in the Green’s operator. We have chosen a strength of the periodic potential that is 
weaker than in present experimental setups. This improves the quality of our approximations, but underestimates 
the size of the spectator state shift in those setups. The spectator state shift grows with the squared amplitude of 
the periodic potential. For example, in setup (A) with a length of 24 cm and a periodic potential strength chosen for 
maximum transition probability for the 1 —> 3 transition, the frequency shifts are a factor of 4 larger. If the length 
was only 10 cm, and again the interaction strength is chosen for maximum transition probability, the frequency shifts 
are about a factor 25 larger, that is, the correction of the resonance frequency of the 1 —> 3 transitions is about 2%. 
Furthermore, we note that for short setups and therefore for a large amplitude excitation, the correction terms in Eq. 
may not be small enough, and the formalism cannot be used. 

The periodic potential in setup (C) appears to be different, as was shown in Refs. [T21 [22] ; we have Wibr = 
am n g sin uit ■ 


■ iahui cos cot-M- with a vibration amplitude a. The first summand does not depend on position, and like 
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TABLE I. Frequency shift for maximum disappearance probability from initial state no with a periodic potential in resonance 
with no —>• n\. 


Transition 

Unperturbed transition frequency (E ni — E no )/2nh 

5 ma .x/27r from Eq. (541 

5max/27r from simulation 

1 -)• 2 

254.564 Hz 

0.28 Hz 

0.24 Hz 

1 -> 3 

462.979 Hz 

-0.03 Hz 

-0.04 Hz 

1 —> 4 

647.177 Hz 

0.16 Hz 

0.17Hz 

1 ->• 5 

815.557Hz 

-0.29 Hz 

-0.29 Hz 

1 6 

972.459 Hz 

-0.018 Hz 

-0.019Hz 

1 -»• 7 

1120.491Hz 

0.45 Hz 

0.4 Hz 


a uniform magnetic field periodic in time discussed above Eq. (22 1 , 


it does not have a visible effect. 
d 


The second 


summand can be discussed with the method described above, using U v ib r = iahuj Table |II| compares the analytical 
result with a numerical simulation, again with a vibration strength that would cause maximum transition probability 
after 0.1s. The analytical calculation is analogues to the one presented above, and uses matrix elements of the type 
(k\J^\l) that have been derived in [ 22 ]. 


TABLE II. Frequency shift for maximum disappearance probability from disappearance probability from initial state no with 
a vibration in resonance with no —> ni. 


Transition 

Unperturbed transition frequency (E ni — E no )/2nh 

5max/27T from Eq. (541 

5max/27r from simulation 

1 -)• 2 

254.564 Hz 

0.14 Hz 

0.15 Hz 

1 -)• 3 

462.979 Hz 

-0.08 Hz 

-0.09 Hz 

1 —> 4 

647.177 Hz 

0.13Hz 

0.12 Hz 

1 ->• 5 

815.557Hz 

-0.31 Hz 

-0.29 Hz 

1 6 

972.459 Hz 

-0.04 Hz 

-0.03 Hz 

1 ->■ 7 

1120.491Hz 

0.44 Hz 

0.44 Hz 


The formally similar effect in atomic physics, called AC Stark shift, has been discussed in a general way in Ref. 
[24] . In atoms, many transition matrix elements vanish due to selection rules, and the nonzero ones are very much 
smaller than differences between energy levels. Usually, large effects are expected only for spectator states that are 
separated from one of the states in the resonant transition by an amount that is very close to hu. In gravitational 
resonance spectroscopy, many spectator states contribute to the shift, and the largest contribution is often given from 
a spectator state that is close to one of the resonant states, although it is not close to be in resonance. 


V. INTERFERENCE SHIFT 

This section discusses setups (A) and (B), where the initial state is prepared by a step. The assumption that the 
initial quantum state of the incoming neutrons is state 3 is an idealization. We can approximate this situation by 
starting with incoming neutrons to be in the ground state before the step that separates region (1.) and (2.). The 
wave function of the incoming neutrons in region (1.) is Here, we have ipi n {z) oc Ai ((2 — Z\ — A;z s tep)/zo) for 

z > Az s tep, where Az step is the height of the step. Again, the one-dimensional description is justified for step heights 
in the order of tens of micrometers, as the energy differences between all relevant quantum states are so small that 
the change in the horizontal velocity of the neutron during a state transition can be neglected. If the step is taken to 
be ideally sharp, we can use the sudden approximation to compute the coefficients in eq. (j3j) just after the step, at 
the beginning of region ( 2 .): 


OO 

a m (after step) := a m ^(after step) = a m- j.(after step) = J tp* n (z)ip m (z)dz (55) 

Azstep 

Fig. 0 shows the coefficients a m (after step) for different step heights around its planned value of Az step = 21 pm. 
The step height is chosen such that the coefficient 03 (after step) is high compared to the other coefficients. Never¬ 
theless, uncertainties in the determination of the step height, and the step shape, and admixtures to the incoming 
neutron state will cause uncertainties in the coefficients of the quantum states which are discussed in this section. For 
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FIG. 4. Initial coefficients a m (after step). These coefficients are independent of the spin state, and therefore the subscript 
denoting the spin is suppressed. 


Az s tep = 21 fim, we have ai(after step) /a ^(after step) ~ 0.05. We have admixtures of quantum states other than 1 
and 3, however, for a resonance frequency close to the lf>3 transition, they are less important. 

Neutrons traversing the distance from the step to the transition region acquire some phase that depends on the 
state number m. This phase does not change the energy levels or the resonance frequency. However, it does changes 
the time needed for maximum transition, and therefore it might look like a resonance peak shift, depending on the 
setup. In setup (A), the phase difference between state 3 and ground state, a = arg(ct3(0)) — arg(ai(0)), is difficult to 
determine from the experiment geometry, and it might have to be treated as a fit parameter. The interference shift 
cannot be observed in setup (B), as here the phase of the magnetic field in the transition region changes with time. 


The filter in region (2.) will be set such that the rate of accepted neutrons N\(v x 
Using eq. (241, we have a count rate N\ (v x ) after the filter that is given by 


is proportional to |ai(t = L/v x )\ 


Ni(v x ) oc |ui 3 1 |a 3 ( 


23? [% 3 una3(0)ai(0)] 


M 2 M0)| 2 


(56) 


The elements of the matrix iq m are taken at t = L/v x . This is the time the neutrons take to fly through the rotating 
field in the transition region. The first term describes how neutrons that are initially in a pure state 3 are appearing 
in state 1 at the end of the transition region. The last term describes how neutrons that are in a pure state 1 at 
the beginning of the transition region are disappearing from this state. Its contribution is reduced by |ai(0)/a3(0)| 2 , 
and therefore small. The second term is an interference term. It disappears at resonance, and it is at maximum 
for |u; — W 3 3± | = |r^ 3 i±|. Its size is reduced by |ai(0)/a3(0)| relative to the first term. At that level, the resonance 
peaks become asymmetric. This is in addition to the trivial (and small) asymmetry due to the fact that the length of 
the transition region L can be optimized to have maximal state transfer for a certain velocity component v x in the 
incoming neutron beam, but not for all of them. Using eqs. ( |24[ ) and 5uj(v x ) = 2irv x /\ — 3 we can predict the 

count rates at given horizontal velocity component to be 


n§ x ± |a 3 (0)| 2 + 2n 31 , ± Suj{v x ) |a 3 (0)| K(0)| cosa + 6u 2 (v x ) |ai(0)| 2 
Ai(d x ) oc - - - ^77 — , -sm 


Suj 2 (v x ) + Slg 


31 ,± 




+ 


^ 31 ,± |a 3 ( 0 )| |ar (0)| sin a 


n 2 


yJS uj 2 (v x ) ■ 
a i(0)| 2 cos 2 (sj(5 uj 2 (v x ) + d 


sin (yJ5 uj 2 (v x ) + 


2v x 


(57) 


Fig. [5] shows the precision with which the model presented here allows to determine the position of the peaks. The 
agreement with the full simulation is well below the peak width for the peaks for which the length of the transition 
region is optimized for maximum transition. For the small peaks, the combination of interference effects and the 
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3 4 5 6 7 


v x [m/s] 


FIG. 5. Expected yield for different velocity components v x of the incoming neutron beam.The arrows in parenthesis denote if 
the reso nan ce is for the "spin-up" (']') or the "spin-down" (1 ) co mponent. The bullets are obtained in a numerical simulation 
of eqs. (C2 I. The gray lines show the analytical model (eq. ( |57| )). The arrows indicate the velocity components for which the 
transitions occur that are indicated above the arrows. The v x distribution has been neglected in this figure. It is expected to 
emphasize the higher values for v x . 


non-perfect reproduction of the Rabi frequency by the model is a limitation. For "spin-up" neutrons, the 1 o 2 
transition is not well separated from the 2 0 3, and our two-state model cannot capture the dynamics correctly. The 
transition 1 f> 4 is a case where the interference shift is large. 

The maximum of the resonance peak, v XyPea k, is shifted from its unperturbed value Aw( 3 ± /27 t due to the interference 
term. In first approximation, this shift is given by 


Av, 


x,peak 


V x,Peak (l- fl+ ^^COt^^) 

\ a 3 (0)/ \ ^13,± ^ ^ / 


CLl ( 0 ) ^ 31 ,± 

a 3 (0) o/ 13)± 


cos a 


Ql(0) ( ^31, 
fl 3(0) \w' 13t 


± Ln 


fcot 2 - 1 

\ ^13,± A 


^ sin a 


The interference shift is small for full transitions, where we have 

Av x ^p e ak Oi(0) 

eak ^3(0) ^13,± 


cos a 


(58) 


(59) 


Irrespective of the value for a, the correction is usually below 0.1%. It gets large for a given resonance if the length 
of the transition region L is not chosen such that the transition probability is at maximum, that is, if the argument 
of the cotangents is not close to an half-integral multiple of 7r. In Fig. [5] the gray arrow show the position of the 
maximum after taking into account the interference shift. The difference to the expectation without interference (the 
black arrow) is small, and most pronounced at the lf)4 transition, where the transition probability is small since 
the length of the system L is not optimized for that transition. 


VI. SUMMARY 

This paper discusses systematic frequency shifts in the spectroscopic measurement of the energies of gravitationally 
bound quantum states of ultra-cold neutrons. The largest one, the Stern-Gerlach shift, is associated with the dynamics 
of the spin; it is specific to the setup presented here, in which resonant transitions between quantum states are induced 
with an oscillating magnetic field gradient in flow-through mode. The Stern-Gerlach shift is spin-dependent, and 
disappears by taking a suitable average over the spin state. The spectator state shift is present in any scheme of 
gravitational resonance spectroscopy. The interference shift stems from an imperfect preparation of the initial state 
through a step, and is averaged out in many schemes. These frequency shifts need to be discussed and taken into 
account in a gravitational resonance spectroscopy if precision measurements are desired. 
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Appendix A: Magnetic field of parallel rectangular wires 


We will now compute the magnetic field of an infinite amount of parallel rectangular wires with regular spacing 
and current pattern. This is the set of wires in the transition region in setups (A) and (B), except for that the real 
wires have finite length, are interconnected, and we have a finite amount of them. We will solve the problem in three 
steps: First we give an expression for an infinite number of parallel thin wires, then we extend the solution to wires 
with rectangular cross section, and finally we superimpose wire sets with different currents to take into account the 
pattern of currents we have in our setup. 

The magnetic field of a single wire with current /, parallel to the y axis, positioned at (xq,Zq), for a point (a :,z) 
outside the wire, is 


B x {x,z) = 


2tt ( x - x 0 f + (z - z 0 ) 2 


(Al) 


For an infinite amount of wires at distance A, we add the magnetic field of each one. Assuming the n th wire is at 
x = xq + n\, we get 




Z~ Zq 


2n ( X - Oo + ™A )) 2 +(Z~ Z 0 f 


(A2) 


To understand better the x and z dependence of this function, let us calculate the Fourier transform with respect to 
the x variable: 


B x (k,z) = -J= J B x (x,z)e lkx dx . 

— OO 

Taking into account that for z < Zo (the case we are interested in), 

y e lkx dx = -7re too e- |fe(2 “ 2o)l , 


(A3) 


/ 


_ Z - Zq 

(x - X 0 f + (z- Z 0 ) 2 


one obtains 


B x (k,z) = —^Le~ l fe < 2 - 2 °)l 

M 7 2\^2n 


e ikx 0 +ikna _ c ~ \ k (z-z 0 ) \ ^ikxg S — 


n =—oo 


27r 


n- 


We used here the usual expression for the Dirac comb (T = 27t/A): 


1 OO 

^ E 


Q i27r nt/T _ 


= £ 8(t-nT) . 


Therefore, for the magnetic field, we obtain 

OO 

B x (x,z) = --2= J B x (k, z)e~ lkx dk 

— OO 

and 


B x (x, z) = - 


Mo I 
2A 


E 


^2Kn(x-xo)/\ p -2ir\n(z-z 0 )\/\ = _ MO 7 ( 1 , ^ r -2irn\z-z 0 \ / A pnn 2™{x - Xp) 


COS 


(A4) 


(A5) 


(A6) 


(A7) 


(AS) 
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This field represents a series of periodic terms where higher harmonics are suppressed by a factor of e 2n ( z z o)/ x : 


B 


.(*, c) = - A (i + e- 2 'l-»l/* cos AifA + cos AiyA + .. ) 

A y2 A A J 


(A9) 


We can rewrite this as 


D / N _ Mo-7 / 1 
B x [x,z) — + 2_ y 


=s -2im|z-2 0 |/A cog 2?r ( X X o) 


A 


_ 7*0-7 l _|_ g —2-7rn|z —z 0 |/Ag27rm(x—x 0 )/A _|_ g-27rn| 


2 :— 2:0 |/Ap— 2nin(x— xq) / X 


Mo I 
2A 
Mo-7 


1 + 


n—1 

g—27rn|2: —20 |/^027ri(x— xq)/A 


n—1 

e — 27r|z—2: 0 |/A e -2-7ri(ai—a? 0 )/A 


^ 0— 27 r| 2 :—2o|/A027ri(ai—aio)/A ^ 0— 2 n\z— 2o|/A0— 2 / ki{x—Xq)/\ 

sinh27r(z — z 0 )/X 


2A cosh27r(2 — zq)/X — cos2tt(x — Xq)/X 
An analogous expression can be obtained for the ^-component of the field: 

,(x, z) = - A ('e-“l>-»l A sin AiyA + e —1~.1A si „ 4 A - A + .. 
Ay A A 

l-iol sin 27r(a; — Xq)/X 


(A10) 


B, 


2A cosh27r \z — Zq \/A — cos27r(:r — Xq)/X (^--^1) 

The next step is to replace the infinitely thin wires at x = Xq + nX and z = Zq by a wire with rectangular cross 
section of width Ax = 1 mm and height A;z = 1 mm, centered at the same place. For z < Zo — Az/2, we get: 

xq-\-Ax /2 zq-\-Az /2 


B x {x,z) = - 


XAxAz V 2 


Xq — Ax/2 2q-Az/2 


/*0 J (l + g—2 W |z—z'|/A cos 2 *( X ~ X ') 

A 

e-l-s'lAcos^Wfl. 

A 


dz'dx' 


(A12) 


We evaluate this expression, to 
B 


(z, z) = - A 7i + sinh Af cos AiyA si „ A£ 

v 7 A V 2 7r 2 Aa;A0 AAA 

A 2 _ 4 _i i _, n i/x . , 27rAz 47r(a; — x 0 ) . 2nAx 

+ \o o —e 4,r|z Zol/A sinh —-— cos —A--— sin —-- 

(2-7r) 2 ArrAz A A A 

In the same way, we get 

B ,(*, z) = - 2d ( sinh Af si „ 2 A ~ A si „ ^ 

A n z AxAz AAA 

A 2 4 _i_ - i/\ . , 27 tAz . 4:n(x — Xn) . 27rAa: 

.-tj-t—— e -4 I °I /A sinh —— sin —A--A s i n — + 

(2tt) 2 AxAz AAA 


(A13) 


(A14) 


The last step is to add the magnetic field of infinite parallel wires with current I = I 0 a,t x 0 = ^offset = A/4 to the 
magnetic field of infinite parallel wires with current I = Iq/^/2. at Xq = a-’offset + A/ 8 , and so on up to the the magnetic 
field of infinite parallel wires with current / = I$j\f2 at xo = x 0 ff S et + 7A/8. For the choice of currents, the terms that 
rotate with spatial period A/ 4 - 7 T and A/ 67 r cancel. Neglecting the terms that rotate even faster, we obtain: 

= -A° Ab e -s.h-s„|/> sinh Af sin 


A ir 2 AxAz 


A 


A 


27r(a: — A/4) | 1 2tt(x — 3\/8) 1 27r(a: — 5A/8) 

A y/2 X y/2 A 

27r(a:-3A/4) 1 27r(a; - 7A/8) | 1 _ 2?r(a; - 9A/8) 

A y/2 x y/2 A 

= - — ■ ^A_ e -s. | .-s.lAsi„h Ai sin A£ s i„ A? . 

A 7r 2 Aa;A0 AAA 


(A15) 
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For the z component, we get: 

B z (x,z) = - 


„ i /\ , 7tA z it Ax 

e 2 ' *°l' A sinh—-—sm- 


Polo A 2 

A tt 2 AxAz A A 

. 2tt(x — A/4) 1 . 2tt(x — 3A/8) 1 . 2tt(x-5X/8) 

sm---- H- -= sm--- -= sm--- 

A ^2 A ^2 A 

2n(x - 3A/4) 1 ^ 27r(x - 7A/8) , 1 ^ 2n{x - 9A/8)\ 

sin r ,— sin - + — sin ~ I 

A y/2 x V2 A J 

= ^ -lA"- 1 "- 1 ' 1 sinh ’Al sin 1+ cos 2= . 

A n 2 AxAz AAA 


(A16) 


From this, we obtain 


B x = B z 

Px = Pz 


Pzz Pxz 


4^ 0 / 0 A 2 

_ — 27r\z — Zq\/X 

A tt 2 AxAz 
8t th 0 I 0 A 2 

A 2 tt 2 AxAz 
167t 2 /j, 0 /o A 2 


7tAz . 7rAx 
sinh —-— sin —-— f 

A A 

e -2,|7- w |/A sinh EAf sin ^x 


1 niT 


A 


A 

IT A x 


„ 2 a A e- 27r ' z -^l/ A sinh ^ sin — 

A 3 7r 2 Aa;Az A A 


0.6T/m 
w 400 T/r 


(A17) 

(A18) 

(A19) 


The numerical values for the amplitudes of magnetic field and its derivatives are given for Iq = 5 A, A = 1cm, and for 
a location that is 1.5mm below the wires (z = zq — 1.5mm). In our discussions, we have used the numerical values 
for B X ,B Z ,... as independent inputs, to be able to take into account a measured magnetic field which might differ 
due to imperfections in the geometry. 

We have shown in HD that v x , the x component of the neutron velocity, can be taken as constant for each neutron 
for the measurements planned. Therefore, we can write the magnetic field alternatively as a function of time and z 
coordinate using the translation x = v x ■ t. The rotation frequency of the field, as seen by the neutron, is then given 
as w = v x /X. 

For this field, the expansion coefficients used in the text in eq. ([7]) can be further expanded in B x /B q and B z /Bq, 
and are given by: 


a o — B 0 \ 1 + 


4 Bl 


~3Bt + 4B 2 B 2 


64Bq 


ai Bz ( 1 8B 2 + 


Bl 3 Bt - ABlB 2 


64B o 


+ ... 


~ (1.5 + 0.167 + 0.005 + ...) mT , 
(1.0- 0.056-0.003 + ...) mT , 


Bt 


a.2 — — 


Bl 


4B n 16Bn 


+-(-0.167 + 0.018 + ...) mT , 


(A20) 


R PxB x -Zp x B 3 x + 2p x B x Bl + 2/3 z B 2 x B z , l5p x B 5 x - 12p z BtB z - 24/3 x B 3 Bl + l&p z B 2 x B 3 z + 8p x B x Bt 

P o = -r 


2B 0 16 B 0 d 

H-(0.2 + 0.0111 + 0.0019 + ...) T/m 


128B 0 5 


Pi — Pz + 


-p z B 2 x - 2 p x B x B z , 3 p z B$ + 12p x B x B z - l2p z B 2 x B 2 z - 8 p x B x B 


SB q 2 


64B 0 4 


~ (0.6- 0.1 - 0.0093 + ...) T/m , 

ft = PXA + /+} ~ (-0.2 + 0.044 + ...) T/m 

Z-DQ 4 


(A21) 


Appendix B: Equivalence of solutions to Schroedinger equation in adiabatic approximation to the original one 

The Schroedinger equation for a oscillating magnetic interaction potential is given by 

= H 0 \^(z,t)) - p n -B(t)\ip(t)) =ih~\4>(t)) . (Bl) 

We had asserted in the main text that this equation is equivalent to the one that we used as eq. ([5]), that is, in 
adiabatic approximation. In this section we want to derive a condition for that, valid for the conditions of setup (A). 
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Let us define the 
the magnetic field. 


spinor |+) (|—)) as one that describes a neutron spin that is always aligned (anti-aligned) with 
That means that this spinor is an eigen-vector of the magnetic interaction Hamiltonian 


- /2„ • B(z, t ) |±) = e±(z, t) |±) 


(B2) 


Here, z and t play a role as parameters. For B(z,t ) = B x (z,t)x + B~(z,t)z, one easily gets e±{z,t) = ±fj, n B(z,t). 
The spinor has the components 


1 +) 


b> 


B-B z 

y/ B 2 +(B-B z ) 2 
_ Hr _ 

y/BZ + iB-B ,) 2 
B+B- \ 

y/Bl + {B+B z ) 2 

_Hr_ 

y/Bl + (B+B z ) 2 ) 


With these spinors, the general wave function is now 



(B3) 

(B4) 


bW) = <mm) i +)+<t>-{z,t) |-) 


(B5) 


Taking into account that the spinors |+), |—) are orthogonal to each other at any moment t and any position z, we 
get the corresponding system for the new wave function components 

d d 

Horf>+(z, t) + t)4> + (z, t) = ih—<j> + {z,t) + ituj)-{z,t) H^l+) > 

d d 

H 0 <j>-{z,t)-HnB(z,t)<j)-(z,t) = ih— <j>-{z,t)-ih{+\ — |-) . (B6) 

Decoupling of amplitudes that corresponding to different spin projections on the field axis (the adiabatic approxima¬ 
tion) takes place if we neglect the last term in each equation, that is, their coupling. This approximation is valid 
if 


*h 5 w 


< \HnB\ 


(B7) 


The spinors |+), |—) are periodic functions of t, and their periodicity is chosen to match the transition frequency, e.g. 

Wi3.±: 


h 


<b' + > 


1 dB z 
2ao dt 


B z nv x 


(B8) 


Here, ag is the average value of the magnetic field magnitude \B\, as defined in Eq. Q. 
frequency with which the magnetic field changes is small compared to the Larmor frequency 


Eq. (B7l 
2 Li n a 0 fh: 


holds if the 


B z nv x ^ 2 ii n a 0 
a 0 X H 


(B9) 


For this condition to hold for all neutron velocities up to v x = 7 m/s, we need a.$ y/\B z hTTV x /2\yL n ). This is verified 
in a numerical simulation of the full Schroedinger equation (B1 ) including both spin states, demanding Bq > 1.5 mT. 


Appendix C: Numerical simulation of transition region in setup (A): 


We have developed a full numerical simulation of the Schroedinger equation with spin (B11. We use the ansatz: 

OO 

I ip(t)) = ^2 |m f) + a m i(t)e~^ Ernt \m i) . 


(Cl) 


m =1 


Here, |m f) is the m th solution of eq. (B11 without magnetic field to a spin pointing into the +z direction ("spin-up") 
to energy E mi and \m 4-) is the m th solution of eq. (B1) without magnetic field to a spin pointing into the -z direction 
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("spin-down"). We define ui m i = 


(Ei — E m )/h. The coefficients a m f(t) and a m i(t ) have to fulfill the coupled equations 


d&m'X 

dt 


da m _. 

dt 


OO . OO 


Z=1 

. OO 

/i n * B 

; t)-^Ea 4 e-^(TOt 

1 1=1 

OO 

/in ‘ B 

j i ^2 a n e ~ lulmlt (™4 

/in ‘ B 

u) - ^ E a n e ~ lUmlt ( m 4- 

“/in ’ B 


(C2) 


For the computation of the matrix elements, we need to write —fl n ■ B = /r n (? • B , where a x , a y , and < 7 - are the 
Pauli matrices. For the piece of the magnetic field linear or quadratic in z, we have used: 


(m t | z\l t) = (m| \z\l I) 
and (to f | z 2 \l f) = (to 4- \z 2 \l 4-) 


(2/3 )zm for to = l 

2(—1 ) m ~ l+1 Zq/( z m ~ zi) 2 otherwise 

(8/15)z^ for to = l 

24(—l) m-Z+1 2 g/(z m - zi) 4 otherwise 


(C3) 

(C4) 


These identities have been derived in f25j, |2B] ■ The sign of the results for to 7 ^ l depends on the sign convention for 
the normalization constants in ip m (z) and (/> i(z ). Throughout this paper, we take the normalization constant to be 
real and positive for each state. 

The numerical simulations in Figs. [2] and [3] are obtained by solving these equations. Note that in both figures, 
we compare the numerical simulation for the population of the ground state to H 0l |ai+| + |aii | , with the model 

prediction for the ground state of Ffo,±, + I®/!”- Due the smallness of (1 f |to±), (1 4- |to±) for m > 1, the 

different meaning of these populations is not visible. 


[1] V.V. Nesvizhevsky, H.G. Bonier, A.K. Petukhov, H. Abele, S. Baefiler, F.J. Ruefi, T. Stoferle, A. Westphal, A.M. Gagarski, 
G.A. Petrov, A.V. Strelkov, Nature 415, 297 (2002). 

[2] V.V. Nesvizhevsky, H.G. Borner, A.M. Gagarski, A.K. Petukhov, G.A. Petrov, H. Abele, S. Baefiler, G. Divcovic, F.J. 
Ruefi, T. Stoferle, A. Westphal, A.V. Strelkov, K.V. Protasov, A.Yu. Voronin, Phys. Rev. D 67, 102002 (2003). 

[3] V.V. Nesvizhevsky, A.K. Petukhov, H.G. Borner, T.A. Baranova, A.M. Gagarski, G.A. Petrov, K.V. Protasov, A.Yu. 
Voronin, S. Baefiler, H. Abele, A. Westphal, L. Lucovac, Europ. Phys. Journ. C 40, 479 (2005). 

[4] S. Baefiler, Journ. Phys. G 36, 10 (2009). 

[5] V.V. Nesvizhevsky, Phys.-Uspekhi 53, 645 (2010). 

[6] V. Nesvizhevsky, K. Protasov, in: D. Moore (Ed.), Trends in Quantum Gravity Research, Nova Science Publishers, New 
York, 2005, pp. 65-107. 

[7] T. Jenke, P. Geltenbort, H. Lemmel, H. Abele, Nat. Phys. 7, 468 (2011). 

[8] T. Jenke, G. Cronenberg, J. Burgdorfer, L.A. Chizhova, P. Geltenbort, A.N. Ivanov, T. Lauer, T. Lins, S. Rotter, H. Saul, 
U. Schmidt, H. Abele, Phys. Rev. Lett. 112, 151105 (2014). 

[9] M. Kreuz, V.V. Nesvizhevsky, P. Schmidt-Wcllenburg, T. Soldner, M. Thomas, H.G. Borner, F. Naraghi, G. Pignol, K.V. 
Protasov, D. Rebreyend, F. Vezzu, D. Forest, P. Ganau, J.M. Mackowski, C. Michel, J.L. Montorio, N. Morgado, L. Pinard, 
A. Renrillieux, S. Baefiler, A.M. Gagarski, L.A. Grigorieva, A.M. Kuzmina, A.E. Meyerovich, L.P. Mezhov-Deglin, G.A. 
Petrov, A.V. Strelkov, A.Yu. Voronin, Nucl. Inst. Meth A 611, 326 (2009). 

[10] G. Pignol, Thesis, Universite Joseph Fourier, Grenoble (2009). 

[11] S. Baefiler, M. Beau, M. Kreuz, V. N. Kurlov, V.V. Nesvizhevsky, G. Pignol, K.V. Protasov, F. Vezzu, A.Yu. Voronin, 
Compt. Rend. Phys. 12, 707 (2011). 

[12] V.V. Nesvizhevsky, A.K. Petukhov, H.G. Borner, T. Soldner, P. Schmidt-Wellenburg, M. Kreuz, G. Pignol, K.V. Protasov, 
D. Rebreyend, F. Vezzu, D. Forest, P. Ganau, J.M. Mackowski, C. Michel, J.L. Montorio, N. Morgado, L. Pinard, A. 
Renrillieux, A.M. Gagarski, G.A. Petrov, A.M. Kusmina, A.V. Strelkov, H. Abele, S. Baefiler, A.Yu. Voronin, Proceedings 
of "15th International Seminar on Interaction of Neutrons with Nuclei", Dubna, Russia (2007) (arXiv:0708.2541). 

[13] I. Antoniadis, S. Baefiler, M. Buchner, V.V. Fedorov, S. Hoedl, A. Lambrecht, V.V. Nesvizhevsky, G. Pignol, K.V. Protasov, 
S. Reynaud, Yu. Sobolev, Compt. Rend. Phys. 12, 755 (2011). 

[14] P. Brax, G. Pignol, Phys. Rev. Lett. 107, 2011 (111301). 

[15] A. Peters, K. Chung, S. Chu, Nature 400, 849 (1999). 

[16] E.G. Adelberger, J.H. Gundlach, B.R. Heckel, S. Hoedl, S. Schlamminger, Prog. Part. Nucl. Phys. 62, 102 (2009). 

[17] V.V. Nesvizhevsky, H. Borner, A.M. Gagarski, G.A. Petrov, A.K. Petukhov, H. Abele, S. Baefiler, T. Stoferle, S.M. 
Soloviev, Nucl. Inst. Meth A 440, 754 (2000). 

[18] G. Pignol, S. Baefiler, V.V. Nesvizhevsky, K. Protasov, D. Rebreyend, A. Voronin, Adv. High En. Phys. 2014, 628125 
(2014). 

[19] H. Abele, T. Jenke, H. Leeb, J. Schmiedmayer, Phys. Rev. D 81, 065019 (2010). 













20 


[20] J.H. Shirley, Journ. of Appl. Phys. 34, 783 (1963). 

[21] F. Bloch, A. Siegert, Phys. Rev. 57, 522 (1940). 

[22] A. E. Meyerovich, I. V. Ponomarev, Phys. Rev. B 65, 155413 (2002). 

[23] A.Yu. Voronin, H. Abele, S. Baefiler, V.V. Nesvizhevsky, A. K. Petukhov, K.V. Protasov, A. Westphal, Phys. Rev. D 73, 
044029 (2006). 

[24] J.H. Shirley, Phys. Rev. 138, 979 (1965). 

[25] D.M. Goodmanson, Am. J. Phys. 68 , 866 (2000). 

[26] R.J. Gordon, J. Chem. Phys. 51, 14 (1969). 



